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Using the (dynamical) density-matrix renormalization group method, we study the low-energy 
physics of three-leg antiferromagnetic Heisenberg model where the periodic boundary conditions 
are applied in the rung direction. We confirm that the spin excitations are always gapped as long 
as the intra-ring couplings form a regular triangle. From precise finite-size-scaling analyses of the 
spin gap and dimerization order parameter, we also find that the spin gap is collapsed by very 
small asymmetric modulation of the intra-ring couplings. Moreover, the dynamical spin structure 
factors on the intra- and inter-leg correlations are calculated. It is demonstrated that the low-lying 
structure of the inter-leg spectra is particularly affected by the asymmetric modulation. 

PACS numbers: 75.10.Jm, 75.30.Kz, 75.40.Gb, 75.40.Mg 
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I. INTRODUCTION 

For many years spin ladder systems have attracted 
much attention not only due to the existence of a vari- 



ety of related materials, e.g., Sr„_iCu„+i02n (Refs.[l|2), 
La4+4nCu8+2TiOi4-i-8« (Rcf. H), and CaV205 (Ref. H), 
etc., but also as intermediates between onc-dimcnsional 
(ID) and two-dimensional quantum spin physics. It has 
been confirmed both experimentaUy and theoreticaUy 
that spin-i ladders are gapful for an even number of legs 
and whereas gapless for an odd number of legs when the 
open boundary conditions (OBC) are apphed in the rung 
direction (e.g., as a review, see Ref. M)- On the other 
hand, if the periodic boundary conditions (PBC) are ap- 
plied in the rung direction (referred as a spin tube) for 
odd-leg ladders, the spin states are drastically changed by 
associating with the occurrence of frustration; it is known 
that the system is spontaneously dimerized to remove the 
frustration and all the spin excitations are gapped.^^ 

At present, there are two experimental candidates for 
odd-leg spin tubes. One of them is vanadium oxide 
Na2V307, which may be regarded as a nine-leg Heisen- 
berg spin tube systcmi^ The ^'^Na NMR response, dc- 
and ac-magnetic susceptibilities, and the specific heat 
measurements^ reveal that above 100 K the system is 
considered as paramagnetic; below 100 K most of the lo- 
calized V magnetic moments {S = ^) form a collection of 
spin-singlet dimers with gaps A ~ — 350 K and the re- 
maining small fraction of them forms spin-triplet bound 
states with gaps A '--^ — 15 K; and the degeneracy of 
the triplet ground states is lifted by a phase transition at 
0.086 K. The low-energy model Hamiltonian of Na2V307 
has been proposed by some theoretical groups. We seem 
to have reached a consensus that the intra-ring exchange 
interactions are antiferromagnetic. However, inter-ring 
ones are still controversial: the ab initio microscopic anal- 
ysisiS argued that they are frustrated antiferromagnetic 
and the magnitude is much smaller than the intra-ring 



ones; in contrast, the first-principle calculationsii esti- 
mated them to be ferromagnetic and of the same order 
of magnitude with the intra-ring ones. 

The other experimental candidate is three-leg com- 
pound [(CuCl2tachH)3Cl]Cl2, which is composed of al- 
ternating (CuCl2tachH)3 triangles along the crystallo- 
graphic c axis The effective model has been considered 
to be a linearly coupled triangle spin rings with antifer- 
romagnetic intra-ring couplings and two frustrating an- 
tiferromagnetic inter-ring couplings. The high-field mag- 
netization measurements suggested that all the couplings 
are of the same order of magnitude; in this situation, it 
was numerically confirmed that the effective model has a 
spin-gapped ground stateJ^ 

In this paper, motivated by such developments in the 
field, we study the low-lying excitations of three-leg an- 
tiferromagnetic Heisenberg spin tube. We assume that 
the fundamental low-energy physics of any odd-leg spin 
tube can be essentially epitomized by that of the three- 
leg spin tube. So far, several theoretical researches have 
been reported for the three-leg spin tube system: primar- 
ily, the bosonization study proposed that the three-leg 
spin tube has a spin-gapped ground statci^ It was nu- 
merically confirmed and found that the system is com- 
pletely dimer-ordered with a broken translational sym- 
metryji It has been also suggested that the spin gap is 
suppressed very rapidly with a lattice modulation in the 
rung directionJ^ Additionally, the system in a magnetic 
fi(; yi5,i6,i7,i8 g^j^jj with some kinds of frustrated inter- 
actionsiSiSSii^ have been examined (the latter model is 
applicable to [(CuCl2tachH)3Cl]Cl2). Thus, the ground- 
state properties of the three-leg spin tube are relatively 
well understood. However, only few theoretical studies 
have been devoted to the dynamic a^^i^^ due to the diffi- 
culty of calculating the dynamical quantities. We there- 
fore calculate the dynamical spin structure factor us- 
ing the dynamical density-matrix renormalization group 
(DDMRG) method)^ which has been successfully ap- 
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FIG. 1: Lattice structure of three-leg Heisenberg model. Ex- 
amples of on-leg and on-rung spin-singlet pairs are shown. 

plied to the ID Heisenberg model of latei^ Before the 
dynamical calculation, the spin gap and the dimeriza- 
tion order parameter are investigated to provide a deeper 
insight into the ground state. It allows us to precisely an- 
alyze the dynamical spin structure factor. Based on the 
results of the static and dynamical quantities, we pri- 
marily discuss the effect of the lattice modulation in the 
rung direction on the ground-state and low-lying excited- 
states properties. 

This paper is organized as follows. In Sec. II, we de- 
fine the three-leg antiferromagnetic Heisenberg tube and 
explain the applied methods for the calculations. The 
effective Hamiltonian for the asymmetric case is also de- 
rived. In Sec. HI, we present the calculated results. The 
coupling-strength dependence of the spin gap and the 
influence of the coupling modulation on the low-energy 
physics are discussed. Section V contains summary and 
conclusions. 



II. MODEL AND METHOD 

A. Hamiltonian 

We consider the three-leg antiferromagnetic Heisen- 
berg model, the Hamiltonian of which is given by 

3 

a = l i a(^a') i 

where Sa,i is a spin-i operator at rung i and leg a(= 
1,2,3). J (> 0) is the exchange interaction in the leg 
direction and Jaa' (> 0) is the exchange interaction 
between the legs a and a' [see Fig[l]. When Jaa' = 
const. (Va, a'), we call it a "symmetric case" and set as 
Ji2 = J23 = ^31 = J±', otherwise, a "asymmetric case". 
We take J = 1 as the unit of energy hereafter. 

B. DMRG method 

We employ the DMRG technique which is a powerful 
numerical method for various (quasi) ID quantum sys- 
temsi^ For the calculation of static properties, we use the 
standard DMRG method and the OBC are applied in the 



leg direction. It enables us to calculate ground-state and 
low-lying excited-states energies as well as static quanti- 
ties quite accurately for very large finite-size systems (up 
to ~ 0(1000) sites). We are thus allowed to carry out 
an accurate finite-size-scaling analysis for obtaining the 
energies and quantities in the thermodynamic limit. For 
each calculation, we keep m = 400 to 2400 density-matrix 
eigenstates in the renormalization procedure and extrap- 
olate the calculated quantities to the limit m — s- cx). 
We note that the m-extrapolation is mandatory in the 
present system ([T]) because our DMRG trial state ap- 
proaches slowly to the exact one with increasing m due 
to very strong spin frustration. In this way, the maxi- 
mum truncation error, i.e., the discarded weight, is less 
than 1 X 10"'', while the maximum error in the ground- 
state and low-lying excited states energies is less than 
10^' — 10^^. For all calculations of the static quantities, 
we study the ladders with several kinds of length L = 24 
to 312 and then perform the finite-size-scaling analysis 
based on the system-size dependence of the quantities. 

For the calculation of dynamical properties, we use the 
DDMRG method which is an extension of the standard 
DMRG method and has been developed for calculating 
dynamical correlation functions at zero temperature in 
quantum lattice modelsi^ We now apply the PBC for 
both the leg (x) and rung (y) directions. With the PBC, 
the system size must be restricted practically up to about 
a hundred but the result is numerically exact because the 
spin operators Si. can be precisely defined by 

with momentum q = {2ttzi/ L,2ttz2/3) for integers 
—L/2 < zi < L/2 and Z2 = —1,0,1. The sum runs 
over all sites of the system. Since the exact definition 
of the momentum-dependent operators with the OBC is 
quite difficult, it would be better to choose the PBC for a 
quantitative estimation of the spectrum. In the DDMRG 
calculation, a required CPU time increases rapidly with 
the number of the density-matrix eigenstates so that we 
would like to keep it as few as possible; meanwhile, the 
(D)DMRG approach is based on a variational principle 
so that we have to prepare a 'good trial function' of the 
ground state with the density-matrix eigenstates as much 
as possible. Therefore, we keep m = 1200 to obtain 
true ground state in the first ten DDMRG sweeps and 
keep m = 400 to calculate the spectrum for ladders with 
length L = 24. In this way, the maximum truncation er- 
ror, i.e., the discarded weight, is about 1 x 10"''^, while the 
maximum error in the ground-state and low-lying excited 
states energies is about 10~^. 

C. Effective model 

For the symmetric case, in the strong coupling limit 
J± ^ J the zero-th order approximation is obtained by 
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an isolated three-spin triangle. It can be easily diagonal- 
ized as the higher-energy eigenstates, 

ITTT), Uii), 
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(3) 



with energy 3Jj^/4 and the lower-energy eigenstates, 
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with energy — 3Jj_/4 where uj = exp(27ri/3). From the 
above four degenerate lower-energy states, we can con- 
struct the effective Hamiltonianf^ii 
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where the chirality operators are defined as t+ = 
\L){R\ andr" = \R){L\, respectively. A previous DMRG 
study confirmed that the effective Hamiltonian ^ has a 
spin gap. which is estimated as 0.277J.— 

When the asymmetry is introduced (J12 = J31 7^ 
the higher-energy states (O are still degenerate and the 
eigenstates of the anisotropic three-spin triangle spin sys- 
tem with energy (2Ji2 + J23)/4. While, the four degen- 
erate lower-energy states (|3]) are resolved into the two 
doublets: ones are 



|1) 
|2) 



\rL)-u;\rR) 

\lL)-cj\iR) 
V2 



(6) 



with the eigenenergy — 3J23/4, and the others are 

^L)+Lo\rR) 



|3> = 
|4) = 



V2 

lL)+cj\iR) 
V2 



(7) 



with the eigenenergy — J12 -I- J23/4. Thus, the effective 
Hamiltonian ^ is modified as 



Since the effect of the asymmetry is written as the 
crystal-ficld-type term, we suppose that the the physical 
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FIG. 2: Spin gap Ag as a function of Jx in the symmetric 
case. Inset is an extended figure for Jx < 10. 



situation is symmetrical with respect to J23 = J12 as a 
function of J23 . Note that the splitting of the four lower- 
energy states ([4]) into the doublets ([6]) and ^ is not com- 
plete for nonzero J (see Appendix). In the limit that the 
crystal-field-type term is large, the effective Hamiltonian 
can be reduced to the 5=1/2 Heisenberg Hamiltonian 
J X]j=i ^3 ■ "Sj+ij which has no spin gap. 



III. RESULTS 
A. Spin gap 

1. symmetric case 

First, we study the spin gap in the symmetric case. 
It is known that for all positive J± the system is spon- 
taneously dimerized and all the spin excitations are 
gappedi^ii It means that all spins contribute to the for- 
mation of spin-singlet pairs. If we could pick up a dimer- 
ized pair of triangles including six spins, three spin- 
singlet pairs are found: one of them is formed in either 
one of the three legs and the other four spins form a 
couple of pairs in the two rungs. As a result, there ex- 
ist two types of spin-singlet pairs, namely, on-leg and 
on-rung pairs. Examples of the on-leg and the on-rung 
spin-singlet pairs are shown in Fig. [T] 

It would be very intriguing to see the quantitative de- 
pendence of the spin gap on J±. The spin gap is eval- 
uated by an energy difference between the first triplet 
excited state and the singlet ground state, 

A,{L) = Ei{L) - EoiL), A, = lim A,(L), (9) 

L — >oo 

where En{L) is the n-th eingenenergy [n ^ corresponds 
to the ground state) of the system with length L, i.e., 
i X 3 ladder. Note that the number of system length 
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must be taken as L = 21, with l{> 1) being an integer 
to maintain the total spin of the ground state as S" = 0. 
In Fig. O we show the DMRG resuhs of the spin gap 
as a function of J±. The plotted values are extrapo- 
lated to the thermodynamic limit L ^ oo (for example, 
the extrapolation scheme for = 1 is demonstrated in 
Fig.[31). We obtain = 0.254 in the limit of J± = oo. It 
is rather smaller than a value estimated in Ref. because 
the finite-size-scaling analysis is different. 

Roughly speaking, the spin gap increases proportion- 
ally to J_L in the small J± (< 3) regime and keeps almost 
constant in the large J± 10) regime. This behavior 
can be interpreted in terms of different origin of the spin 
gap for each the J± regime, although the mechanism of 
gap opening is invariant for the entire J± regime. In 
other words, the spin gap is approximately scaled by a 
binding energy of most weakly bounded spin-singlet pair 
in the system; and, it switches between the on-leg and on- 
rung pairs at some value of Jj_(~ 5). A more concrete 
description is givn in the following paragraph. 

For J ^ J_L, we can easily imagine that the on- rung 
spin-singlet pairs must be bounded more solidly than 
the on-leg ones. The spin gap is therefore scaled by the 
binding energy of an on-leg pair, i.e., Ao- oc J. Accord- 
ingly, Act is independent of Jj_ and it is consistent with 
the constant behavior of A^ with J± at J± > 10. On 
the other hand, the situation is somewhat different for 
J± < 0{J): the bound state of the on-leg pairs is ex- 
pected to be more solid than that of the on-rung ones. It 
is because that the system is strongly dimerized with in- 
finitesimally small J^. The dimerization strength devel- 
ops abruptly at J_l = 0^ and increases rather slowly with 
increasing Jj_i^ Thus, the spin gap is essentially scaled 
by the binding energy of an on-rung pair. In addition, we 
may assume that the binding energy of the on-rung pair 
is proportional to Jx in the small Jj^ regime, by anal- 
ogy with that of the two-leg Heisenberg systemi^ Now 
therefore, the spin gap is scaled by Jj^, i.e., Ajr oc Jx, 
which is consistent to a linear behavior of A^ with 
at Jx ^ 3. Note that the derivative dAc/dJ± is very 
small 0.053) due to strong spin frustration among the 
intra-ring spins. Consequently, a crossover between the 
constant region and the proportional A^- region is 
seated not at J± ~ 1 but around Jj_ « 5. The existence 
of this crossover can be also confirmed with studying the 
J_L dependence of the dynamical spin structure factor. It 
will be discussed in Sec. IIII Dl 



2. asymmetric case 

Of particular interest is the evolution of the spin gap 
onto an asymmetric modulation of the intra-ring cou- 
plings. Previously, it has been suggested that the spin 
gap is suppressed rapidly by imposing the asymmetry and 
there exists the Berczinskii-Kostcrlitz-Thoulcss (BKT) 
type transition between gapped and gapless phases at 
a finite asymmetric modulationi^i. In order to take the 




FIG. 3: DMRG results of the spin gap as a function of 1/L 
for several values of J23 with fixed J12 = J31 ~ 1. Inset: 
rescaled spin gap as a function of 1/L" with u — 0.649 and 
0.647 for J23 = 0.99 and 1.01, respectively. 

asymmetric modulation into account, we vary the value 
of J23 from unity with keeping J12 = J31 = 1 in our 
model dll). The system-size dependence of the spin gap 
for several values of J23 is shown in Fig. [3l For the sym- 
metric case, J23 = 1, Aa-{L) can be readily extrapolated 
to 1/L ^ and A^ = 0.052(1) is obtained. For a rela- 
tively large asymmetric case, e.g., J23 = 0.9, we can see 
that the spin gap is obviously extrapolated to zero in the 
thermodynamic limit. 

Let us then consider a small asymmetry by making 
J23 only 1% smaller from the symmetric case, i.e., we 
set J23 ~ 0.99. As seen in Fig. [31 Aa-{L) behaves quite 
similarly to that of the symmetric case for small systems 
(L ^ 50); however, the deviation comes into the open 
around L = 50 and increases rapidly as 1/L decreases. 
Even at a rough estimate, Ao-(L) seems to be extrapo- 
lated to a much smaller value at 1/L ^ than that in the 
symmetric case. For more precise extrapolation, a good 
fitting function ought to be chosen. As shown in the in- 
set of Fig.[3l we find that Acr(L) can be scaled better by 
{v < 1) than by in the case of small asymme- 
try. This rescaling allows us a reasonably performance 
of the finite-size-scaling analysis to the thermodynamic 
limit. The spin gap is thus obtained as zero within the 
margin of error ~ ±1 x 10~^. We also find that the 
data for J23 = 1.01 are quantitatively the same as those 
for J23 = 0.99, as expected from Eq. ([S]). These results 
could suggest that the spin gap vanishes as soon as a 
rather (or infinitesimally) small asymmetry in the intra- 
ring exchange couplings is introduced. This is contrary 
to the BKT-type transition suggested in Ref. Q For 
further support of our statement, we evaluate a dimer- 
ization order parameter, which quantifies the presence or 
absence of long-range dimer order indicating the spin- 
gapped ground state, in the following subsection. 
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We should comment on the finite-size-scahng analy- 
sis. The fitting function L~'^ {v < 1) may seem to be 
unusual: however, it is reasonable if wc take into con- 
sideration the fact that the finite-size system will be al- 
ways dimerized due to the Fricdel oscillation. As a conse- 
quence of the dimerization, the spin gap indeed "opens" 
in any finite-size system even if the spin gap closes in the 
thermodynamic limit. We need to remove this "anoma- 
lous" finite-size effect caused by the Friedel oscillation 
to obtain a correct value in the thermodynamic limit. 
Assuming the dimerization strength due to the Friedel 
oscillation decays as (/3 > 0)^^ the spin gap would 
be scaled as L~20/3 ^j^g analogy of the results for 
the ID spin-Peierls Heisenberg model.— Therefore, we 
can justify the finite-size-scaling analysis with the fitting 
function Ac,(L) c>c L~'^ (y = 2/3/3). Incidentally, this 
finitc-sizc-scaling analysis works only for gapless cases; 
actually, the data for the symmetric case cannot be fit- 
ted with A„(L) oc L-". 



B. Dimerization order parameter 

Next, we evaluate the dimerization order parameter 
which indicates the presence or absence of long-ranged 
dimerized state. When the spin gap opens, the sys- 
tem has to be dimerized along the leg direction, i.e, the 
ground state has to be a spin-Peierls one. Therefore, a 
disappearance of this order parameter in the thermody- 
namic limit corresponds to a collapse of the spin gap. In 
order to know whether the spin gap disappears with small 
asymmetry, we study the dimerization order parameter 
in the vicinity of the symmetric case. 

Because the translational symmetry is broken due to 
the Friedel oscillation under the application of the OBC, 
the dimerized state is directly observable with the DMRG 
method. We are interested in the formation of alternating 
spin-singlet pairs in the leg direction, so that we calculate 
the nearest-neighbor spin-spin correlations. 



S{a, z) = - ( Sa,i ■ Sa,i+l 



(10) 



where (• • •) denotes the ground-state expectation value. 
In Fig. lUa) , wc show an example of the Friedel oscillation 
appearing in the nearest-neighbor spin-spin correlations 
of the system with length L = 24 for a symmetric case 
= 1. (In the symmetric case, the results for all values 
of a are the same.) Generally, the Friedel oscillations in 
the center of the system decay as a function of the system 
length. If the amplitude at the center of the system 

D{a, L) ^ |5(a, L/2) - 5(a, L/2 + 1)| (11) 

persists for arbitrarily long system length, it corresponds 
to a long-range dimerization order which indicates the 
spin-Peierls ground state. We thus define the dimeriza- 
tion order parameter as 
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D{a) = lim D{a,L). 

L— i-oc 



(12) 
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FIG. 4: (a) Friedel oscillation in the nearest-neighbor spin- 
spin correlations of the system with length L = 24 for a 
symmetric case J± — 1. (b) Log- log plots of the amplitude 
D{a, L) as a function of the system length L for J23 = 1, 0.99, 
and 0.9 with fixed J12 = J31 = 1. The dotted lines are the 
fitting functions D[a, L) cx L ' (7 > 0). 



In Fig. mb), we show the log-log plots of the amplitude 
D{a,L) as a function of the system length L for several 
values of J23 with fixed J12 = J31 = 1. For the symmet- 
ric case, J23 = 1, the derivative —d\ogD{a,L)/d\ogL 
appears to diminish gradually with increasing L and 
D{a,L) saturates at a value, i.e., £'(0;) = 0.0276(2) in 
the large L limit. It signifies the long-range dimerization 
order, which is consistent with the existence of finite spin 
gap. For reference, we obtain D{a) = 0.0147(5) in the 
large L limit at Jj^ = 0.5 and possibly this value may be 
proportional to J± in the small J± regime. 

Let us then study how the order parameter is affected 
by the asymmetric modulation. For a relatively large 
asymmetric case, J23 — 0.9, the log-log plot of D{a,L) 
with L represents a straight line, i.e., D{a,L) oc L"^-^'* 
(Va), over all ranges of L or at least L ^ 50 [see the inset 
of Fig. HKb)]. We thus find a power-law decay of the 
dimerization order parameter as a function of L. Since a 
power-low decay with the distance indicates the absence 
of long-range order, the disappearance of the spin gap 
is confirmed. We also see that the amplitude D{1,L) is 



6 




FIG. 5: (a) Spin gap Acr{L) as a function of 1/L for several 
values of J23 with J12 — J31 — 100. The solid lines are the 
polynomial fit. (b) Log-log plots of the amplitude D{a,L) 
as a function of L. The dotted lines are the fitting functions 
D{a, L) oc L 7 is estimated as 0.601 for J23 = 99. 



smaller than the others D{2,L)[= D(3,L)] though their 
decay lengths are the same. It is so because the reduction 
of J23 prevents the spin-singlet pairs from forming on the 
rungs between the leg a = 2 and a = 3; accordingly, 
the formation of the pairs on the leg a = 1 is strongly 
suppressed. Note that the (nearly) perfect straight line 
fit of \ogD{a,L) versus logL could be concerned with 
the (almost) linearly scaled spin gap with 1/L. 

We now turn to the case with a small asymmetric mod- 
ulation. J23 = 0.99. The derivative of — logD(Q!, L) with 
logL decreases with increasing L, as seen in the sym- 
metric case. However, it seems to get at a finite value 
around L ^ 150 and the plots for L > 150 can be fitted 
by a straight line. It again implies a power-law decay 
of the order parameter in the large distance. The fitting 
function is estimated as D{a,L) cx £qj. q,^ 

which the slope is much gentler than that in the case of 
J23 = 0.9. Although the decay of the order parameter 
is very slow, the formation of the long-range dimeriza- 
tion order has to be broken down. This slow decay may 
be the reason why the finite-size-scaling analysis of the 
spin gap is quite difficult (see Fig. [3]). Moreover, it must 
be a good guess that the decay ratio of the dimerization 
order parameter develops continuously from zero at the 
symmetric case. Therefore, we argue that the spin gap 
vanishes immediately an infinitcsimally small asymmetry 
is introduced. 

The collapse of the spin gap with small asymmetry 
can be seen more evidently in the strong coupling regime 
Jj_ ^ J. As an example, we show (a) the spin gap as 
a function of 1 /L and (b) the dimerization order param- 
eter as a function of L for J12 = J31 = 100 in Fig. [S] 
For the symmetric case, J23 = 100, the system has a spin 
gap which is estimated as = 0.253 and the dimeriza- 
tion order parameter converges to D(a) = 0.067 in the 
L — > 00 limit. Let us now modulate the triangle rings 
only by 1%, namely, J23 = 99 (and 101). We obviously 
find that the spin gap is extrapolated to zero in the ther- 
modynamic limit and the dimerization order parameter 
decays as a power law with the system size. 



C. A simple intuition 

It would be important to provide an intuitive under- 
standing of the collapse of the spin gap with a small 
asymmetry. Let us consider an isolated triangle Heisen- 
berg ring. We assume the three coupling constants to be 
K, K', and K, corresponding to J12, J23, and J31 in our 
model, respectively. For K — K', the ground state of 
the ring is four-fold degenerate and the spins are com- 
pletely frustrated. When K ^ K', the degenerate states 
are splitted into two states and the spin gap A opens. If 
K' > K {K > K'), a spin-singlet (spin-triplet) pair with 
A = K' — K {A = K — K') between two sites coupled 
by K'. The spin-singlet state corresponds to the doublet 
ones |1). |2) [Eq.Q]; whereas, the triplet state corre- 
sponds to the doublet ones |3), |4) [Eq,®]. We thus find 
that a bound state is stabilized for K ^ K' . We now 
go back to the three-leg tube. Originally, the sponta- 
neous dimerization is vital to form the bound states in 
the symmetric case. As a results, the spin degrees of free- 
dom arc quenched and the spin frustration is weakened. 
However, if the asymmetry is introduced, the bound state 
is naturally formed in each triangle as mentioned above. 
Consequently, the dimcr order is no longer necessary. By 
the analogy with results of the single triangle, the bind- 
ing energy of the bound state in each triangle is scaled 
by |Ji2 — J2'i\ (or IJ31 — J23I) in the asymmetric case. 
This is consistent with the fact that the effective Hamil- 
tonian ([8]) includes the effect of the asymmetry only as 
the crystal-field-type term. 

D. Spin excitation spectra 

1. symmetric case 

Finally, we study the dynamical spin structure factor 
to investigate the low-energy excitations. For the sym- 
metric case, the spin structure factor is defined as 

S{q,uj) - ^(^o|^%HV^„)(V'n|^||V^o) 

n 

x5{u-E^ + E^), (13) 

where \^n) is the 71-th eigenstate with the eigenenergy 
En- The operator {qy = 0, ±|7r) is the Fourier trans- 
formation of the spin operator Sf at site i and given by 
Eq.(l2]) when the PBC are applied in the leg direction. 
In = where the system consists of three separated 
chains, the ground state of the whole system is the di- 
rect product of the ground state of three chains and the 
dynamical spin structure factors are equivalent to those 
of the ID Heisenbcrg model*^'^'^ thus, the spectra for 
qy ~ and qy = ^TT are equivalent. When a finite Jj_ 
is taken into account, the frustration arises among the 
intra-ring spins. In Fig. [51 we show the DDMRG results 
of S{q,u!) at J_L = 1 for a ladder with length L = 24. A 
pronounced peak is found at g = (tt, |7r). It implies that 
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FIG. 6: Dynamical spin structure factor S{q, u)) at the sym- 
metric case J_i_ — \ for (a) qy — Q and (b) q — |7r. The 
system size is fixed at L = 24 and the broadening ry = 0.1 is 
introduced. Insets: S{q,uj) at = tt for J± = 1 (no symbol), 
2 (squares), 5 (circles), and 10 (crosses). 



the spins approximately form a (nearly) 120° structure 
in each triangle ring and an antiferromagnetic correla- 
tion is dominant along the leg direction. On the other 
hand, the peaks around q = (tt, 0) are fairly suppressed 
as compared with the case of Jx = 0. 

For finite J^, all the spin excitations are gapped. 
The position of the lowest-lying peak at q = (tt, 0) 
\q = (tt, Itt)] essentially corresponds to the energy loss 
to break an on- leg [on-rung] spin-singlet pair. Thus, 
the lowest-excitation energy at either q = (tt, 0) or 
q — (tt, Itt) provides the spin gap A^.. For example, 
since the on-leg spin-singlet pair has larger binding en- 
ergy than the on- rung one at J_l = 1, the lowest- lying 
peak of g = (tt, ^tt) is located at lower frequency than 
that of g = (tt, 0) in Fig. [Sj As a result, the position 
of the largest peak in g = (tt, |7r) stands at the low- 
est spinon excitation whose frequency corresponds to the 
spin gap. The magnitude relation of the binding ener- 
gies between the on-leg and the on-rung pairs switches 
at w 5 as mentioned in Sec. IIII A| so that the mo- 
mentum q giving the lowest excitation is also expected 
to change at w 5. 

Let us therefore investigate the evolution of S{q, ui) at 



(0,- 



and q = (tt, -^tt) with J±. The results are 



shown in the insets of Fig. [B) Both the position and the 
weight of the lowest-lying peak at g = (0, |7r) are not that 
much sensitive to J± , which reflects that the dimerization 
strength is almost saturated above J± = 1; whereas, at 
g* = (tt, Itt) the position of the peak shifts towards higher 
frequencies and the weight goes down with increasing J± . 



We then find that the peak position of g = (tt, ^tt) comes 
up with that of g = (0, |7r) at J± « 5. It means that 
the momentum of the lowest-lying peak is changed from 
g = (tt, Itt) to g = (0, |7r) at J± « 5. In fact, this change 
of the momentum corresponds to the switch in the origin 
of the spin gap from the on-rung pair to the on-leg one. 

We now refer a couple of the other significant features: 
one is the appearance of two bands around lu ~ 1.5 and 
3.0 in the g^ = |7r spectra, which are denoted as = cji 
and UJ2 in Fig. [6l respectively. Those bands are associ- 
ated with internal excitation on each triangle ring. Con- 
sidering an isolated triangle ring, there are three states; 
namely, a spin-| quadruplet with energy | [Eq. ^] and 
two degenerate spin-i doublets with energy — | [Eq. dH)]. 
Hence, when J = the excitation spectra show two flat 
bands at wi = and uj2 = ^J±- Thus, we can eas- 
ily guess the spin structure factor in the strong coupling 
regime J± 3> J: the spectra of the low-energy excitations 
[oj ^ 0{J)] are essentially the same as those of the ID 
spin-Peierls (SP) Heisenberg model^'^'^'^^'^; and sin- 
gle peaks associated with the intra-ring excitation exist 
around lu2 = J + ^J± for all of qx- The other feature is 
the reduction of the apparent width of the continuum in 
the qy = Itt spectra, relative to that of the ID Heisen- 
berg model. In other words, the peaks at the lower edge 
of the spectra are rather the (5-function like. It possibly 
reflects the local excitation of the spinon bound state on 
each dimer, as in the spectra of the ID SP Heisenberg 
modeli^ 



2. asymmetric case 

We then examine the dynamical spin structure factor 
in the asymmetric case. Particularly, we focus on the 
change of the spectra with the asymmetry. Since the 
translation symmetry in the rung direction is no longer 
present for the asymmetric case, we redefine the dynam- 
ical spin structure factor as 

5""'(g,,c.) = ^(^ol^^.^^JV^.X^.I^^'.^J^o) 

xS{lj-E, + Eo) (14) 

with 



\ ptq^x qz 



(15) 



where there are four independent combinations of a and 
a', i.e., (a, a') = (1, 1), (2, 2), (1, 2), and (2, 3). 

We here consider the system with a 10% asymmetric 
modulation, i.e., J12 = J31 = 1 and J23 ~ 0.9. The 
DDMRG resuhs of S°'^{qx,u;) for a ladder with length 
i = 24 are shown in Fig. [T] The spin gap should close 
at g^ = 0, TT and the intra-leg spin structure factors 
S^^{qx,i^) and S'^'^ (q^ , oj) are basically similar to that of 
the ID Heisenberg modeL '^^'^^i'^^ The lower edge of the 
spectra is well- fitted with a function lo = csing^,, and c 
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FIG. 7; Dynamical spin structure factor S'"''(ga;, tj) at the 
asymmetric case J12 = J31 = 1 and J23 = 0.9. The system 
size is fixed at L = 24 and tlie broadening 77 = 0.1 is intro- 
duced. The dashed lines denote the fitting of the lowest-lying 
peaks with lu ~ sin^ii. 



is estimated as 1.9 and 1.7 for S^^[qx,ijj) and S'^'^ {q^ , io) , 
respectively. One of the most noticeable deviation from 
the spectra of the ID Heisenberg model is the rapid de- 
crease of the spectral weights S{qx) with distance from 
qx ~ T^- It is because that the spin fluctuations are sup- 
pressed by the exchange interaction between the legs and 
the antiferromagnetic correlations in the leg direction are 
'longer-ranged'. Since the leg a = 1 is more strongly 
coupled with the neighboring legs than the leg a = 2, 
the spectral weights of S^^{qx,io) at small qx are fewer 
than those of S'^'^ {qx , oj) ■ For example, the similar fea- 
ture has been confirmed in the spin structure factor of 
the two-dimensional Heisenberg model which has an an- 
tiferromagnetically ordered ground statei^ 

Let us turn to the inter-leg spin structure factors 
S^'^{qx,uj) and S'^^{qx,uj). These factors themselves rep- 
resent changes of the spectral features derived from the 
asymmetric modulation, because they must be zero in the 



symmetric case. In S^'^{qx,uj), negative- weighted peaks 
appear around qx = ■ It means that the the antiferro- 
magnetic correlation between the legs a = 1 and a = 2 
increases from the frustrated 120° spin structure for the 
symmetric case. Whereas in S'^^{qx,oj), smaller peaks 
with positive weights appear around qx ~ tt. It signifies 
the emergence of the ferromagnetic correlation between 
the legs a = 2 and a = 3, which is caused by the su- 
perexchange interaction via the leg a = 1. With increas- 
ing the asymmetric modulation, the low-energy physical 
properties of the three-leg ladder with the PBC in the 
rung direction seem to quickly approach to those of a 
three-leg ladder with the OBC^ 



IV. SUMMARY 

We study the low-lying excitations of the three-leg 
antiferromagnetic Heisenberg tube with the (D)DMRG 
method. For the symmetric case, we argue that the 
spin gap is scaled by the binding energy of the on- 
rung spin-singlet pair in the weak-coupling regime ( J_l ^ 
3); whereas, the on- leg spin-singlet pair in the strong- 
coupling regime (J_l ^ 10). We then take an asymmetric 
modulation of the intra-ring exchange couplings into ac- 
count. For small asymmetries, precise finite-size-scaling 
analyses of the spin gap and dimerization order parame- 
ter are carried out. Based on the results, we suggest that 
the spin gap vanishes as soon as an infinitesimally small 
asymmetry is introduced. 

Furthermore, we calculate the dynamical spin struc- 
ture factors. In the symmetric case, the low-energy spec- 
tra are essentially the same as those of the ID SP Heisen- 
berg model. Note, however, that additional peaks asso- 
ciating the intra-ring excitation exist. In the asymmetric 
case, the intra-leg spectra are basically similar to those 
of the ID Heisenberg model. They are hardly affected 
by the asymmetric modulation except the spin-gap clos- 
ing at the band edge. On the other hand, the inter-leg 
spectra are profoundly affected even by a small asym- 
metric modulation. An enhancement of ferromagnetic 
or antiferromagneic correlations between two legs can be 
clearly seen. It means that the low-energy physics of the 
three-leg Heisenberg tube approaches quickly to that of 
the non-tube three-leg ladder with increasing the asym- 
metry. 
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FIG. 8: Ground-state expectation values of the projection 
operators (Pi> (= (P2)) (empty) and {P3) (= {P4)) (filled) 
for J12 = J31 — 1 (squares) and 100 (circles). 

APPENDIX 

In this Appendix, wc consider how effectively the four 
degenerate states dH) are resolved to either of the dou- 



blets (O or ([7l) in the ground state. If only configurations 
containing one of the doublets are included in the asym- 
metric case, the Hilbert space to be considered can be 
much reduced. This reduction is exact only in the limit 
of J — *■ and thus we check the validity of the reduction 
numerically for nonzero J. 



When I J12 — J23I ^ J, the ground state of each tri- 
angle is approximately expressed as a linear combination 
of the states |1) and |2) for J23 > J12; whereas, |3) and 
|4) for J23 < Ji2. However, the four states are equally 
likely 'mixed' around J12 = J23 for nonzero J. We here 
prepare projection operators Pi = (i = 1 • • - 4). In 

Fig. (Ha), the ground-state expectation values of the pro- 
jection operators (Pi) (= (P2)) and (P3) (= (P4)) as a 
function of \Ji2 — J23I are shown for a couple of cases 
J12 = J31 = 1 and 100. We can see that the domi- 
nance of the doublet state exchanges continuously be- 
tween |1),|2) and |3),|4) at |Ji2 - J23I ^ 0{J/2). It is 
confirmed for J12 ~ J31 ~ 100 that either (Pi) or (P3) is 
~ 1/2 and the other is ~ at |Ji2 — J23I ^ J/2, where 
the reduction to the doublet of the Hilbert space in each 
triangle must be applicable. Since we have the relation 
-Pi - -P2 + ^3 + P4 = wV+ + WT", the values (Pi) 
and (P3) are symmetric on the reflection against the line 
J23 — J12 = 0, which is consistent with the crystal-field- 
type term in the effective Hamiltonian ([5]). We also plot 
the same results as a function of J23/J12 in Fig. [D^b). 
For J12 = J31 = 1, the situation seems to be more com- 
plex but the properties are expected to be qualitatively 
similar to those in the strong-coupling limit J12, J23 3> J. 
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